MICROCOPY  RESOLUTION  TEST  CHART 

NATIONAL  BUREAU  OF  STANDARPS-1963-A 


REPORT  DOCUMENTATION  PAGE 


4.  TITl I  (m*  JUMUO 

RESEARCH  ON  ALGORITHMS  FOR  ADAPTIVE 
ARRAY  ANTENNAS 


7i  AUTHOR? 

B.  Widrow 
R.  Gooch 
T.  AuTruon 


>•  PERFORMING  QROAMI Z  ATION  HAMS  ANO  AOONESS 

Stanford  University 

Dept  of  Electrical  Engineering 

Stanford  CA  94305 


It.  CONTROLLING  OFFICE  HAMS  ANO  AOOAUI 

Rome  Air  Development  Center  (DCCD) 
Griffiss  AFB  NY  13441 


I*.  MONITORINO  AGENCY  NAMC  A  AOONKSSflf  aWM  tmm  CammIMm  OlHmm) 


STfttSUTtOM  STATCMCNT  (9t  (Ml  *•»«¥) 


Approved  for  public  release;  distribution  vmlimited. 


17.  DISTRIBUTION  STATCMCNT  (ml  Him  Milmi  «NnC  la  BlaaM  70.  II  A/Imm  Im  CipartJ 


I*.  SUPPLEMENTARY  NOTH 

RADC  Project  Engineer:  John  A.  Grander o  (DCCD) 


!f.  KEY  WORDS  (Claim**  ma  imwmmm  no  if  NCMMrr  MMf  Idaaltty  Or  Nxt  ■ 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


RECIPIENT'S  CATALOO  NUMBER 


S.  TYRE  OR  REPORT  S  RERIOO  COVERED 

Interim  Report 
Feb  81  -  Feb  82 


I.  RERRORSHNO  010.  REPORT  NUMBER 

N/A 


9a  CONTRACT  ON  ON  AN  T  NUMNKNfAJ 


F30602-80-C-0046 


PROGRAM  ELEMENT.  PROJECT.  TASK 
AREA  A  WORK  UNIT  NUMBERS 

61102F 

2305J805 


U.  REPORT  OATS 

November  1982 


SECURITY  CLASS,  (mi  Mm  tmpmtt) 


UNCLASS IF 


Adaptive  Systems 
Adaptive  Filters 
Adaptive  Antennas 


Adaptive  Arrays 
Null  Steering 
HR  Filters 
Lattice  Filters 


«•  AMT  A  ACT  (Cwmtitm*  m  nntn  NMb  if  <mnwy  l+mtiiy  ip  Nnk  m— 

7 This  report  demonstrates  a  method  for  Implementing  a  wlde-band  constrained 
array  processor  using  filters  possessing  both  poles  and  zeros.  The  idea 
of  constrained  array  processing  is  introduced  and  an  optimal  array 
weighting  formula  is  derived.  This  optimal  weighting  is  meant  to  cancel 
m-1  wlde-band  interference  sources  incident  on  an  m-element  array  with  a 
constrained  look  direction.  From  the  optimal  weighting  formula,  it 
becomes  evident  that  the  use  of  filters  possessing  both  poles  and  zeros 
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L  INTRODUCTION 

This  report  demonstrates  a  method  for  implementing  a  wide-band  constrained  array  proces¬ 
sor  using  filters  possessing  both  poles  and  zeros.  In  section  II,  the  idea  of  constrained  array  pro¬ 
cessing  is  introduced  and  an  optimal  array  weighting  formula  is  derived.  This  optimal  weighting 
is  meant  to  cancel  m-1  wideband  interference  sources  incident  on  an  m -element  array  with  a  con¬ 
strained  look  direction.  From  the  optimal  weighting  formula,  it  becomes  evident  that  the  use  of 
filters  possessing  both  poles  and  seros  to  perform  the  frequency-dependent  weighting  has  the 
potential  for  improving  the  wideband  nulling  capability  of  the  array.  A  novel  method  for  adapt¬ 
ing  the  pole-zero  array  filters  in  tap-deby-line  form  using  the  LMS  or  RLS  algorithm  is  presented. 

m 

In  the  remaining  part  of  this  section,  the  issues  of  look-direction  signal  bias  and  stability  of  the 
inverse  filter  are  addressed. 

In  subsequent  sections,  a  method  for  adaptively  implementing  the  pole-zero  array  filters  in 
lattice  form  is  detailed.  It  is  shown  that  the  use  of  an  adaptive  lattice  filter  offers  a  considerable 
improvement  in  the  algorithm's  speed  of  convergence.  Convergence  rate  is  an  important  charac¬ 
teristic  of  an  adaptive  algorithm.  A  rapidly  converging  algorithm  eliminates  its  dependence  on 
initial  conditions  and  improves  its  ability  to  track  time-varying  scenarios.  In  adaptive  array  pro¬ 
cessing,  gradient-descent  algorithms  are  widely  used  because  of  their  simplicity  of  implementa¬ 
tion.  When  the  input  signal  autocorrelation  matrix  R  is  "dose”  to  identity,  the  convergence  rate 
of  a  gradient-descent  algorithm  is  as  fast  as  a  Newton-descent  algorithm.  The  convergence  rate  of 
a  gradient-descent  algorithm  is,  however,  dependent  upon  the  eigenvalue  disparity  of  the  R 
matrix.  On  the  other  hand,  bttice  algorithms  are  known  to  be  rebtiveiy  insensitive  to  eigenvalue 
spread.  This  insensitivity,  along  with  an  otherwise  fast  convergence  capability,  should  reduce  any 
concern  about  the  computational  complexity  of  the  lattice. 

In  section  m,  we  present  the  scalar  lattice  predictor  in  the  context  of  speech  processing 
where  it  was  first  introduced.  We  show  that  the  inverse  of  a  bttice  predictor  can  be  implemented 
in  lattice  form  and  that  its  stability  can  be  easily  checked.  We  also  present  the  joint-process  lat¬ 
tice  used  for  estimating  one  signal  based  on  observations  of  another  signal.  Simulations  of  the 
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algorithm  are  presented  to  show  its  performance  onder  a  variety  of  different  input  eigenvalue 
spreads. 

We  continue  in  section  IV  by  presenting  the  2-channel  lattice  predictor  and  show  how  it  can 
be  used  for  equation  error  pole-iero  filtering  by  embedding  the  primary  and  auxiliary  inputs  into 
a  multi-channel  input.  Simulations  showing  the  2-C  lattice  predktor  performance  when  there  is 
an  eigenvalue  spread  are  presented. 

In  section  V,  we  apply  the  lattice  predictor  to  a  two-element  pole-sero  array  processing 
problem.  The  question  of  forming  an  inverse  filter  to  restore  the  true  output  error  is  solved  via 
the  computation  of  an  equivalent  scalar  lattice  predictor.  To  complete  this  section,  antenna  pro¬ 
cessing  simulations  using  the  pole-sero  generalised  sidelobe  canceler  are  presented. 

In  section  VI  we  present  an  extension  of  the  pole-sero  lattice  predictor  to  multi-element 
antenna  processing.  We  point  out  that  the  computational  complexity  of  such  a  multivariable 
filter  grows  as  the  cube  of  the  dimension  of  the  input  vector  and  directly  with  the  filter  order. 

And  last,  in  section  VII  we  present  our  conclusions  and  suggestions  for  future  research. 
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n.  WIDE- BAND  ADAPTIVE  ARRAY  PROCESSING  USING  POLES  AND  ZEROS 

M  Introduction 

Conventional  adaptive  array  processing  is  accomplished  by  linearly  combining  the  ontpnts 
of  tap  delay  lines  attached  to  each  sensor  of  an  array.  This  type  of  processing  can  be  interpreted 
as  using  an  "all-sero"  filter  at  each  sensor  to  generate  a  frequency-dependent  magnitude  and 
phase  shift  (weighting)  over  the  desired  array  bandwidth.  In  this  section  of  the  report  we  present 
a  new  method  for  performing  the  frequency-dependent  array  weighting  using  filters  possessing 
both  poles  and  seros.  Adaptation  of  these  filters  is  based  on  "equation  error”  rather  than  on  the 
usual  "output  error”.  The  method  is  shown  to  substantially  improve  the  wideband  interference 
nulling  capabilities  of  the  array,  providing  sharper  and  deeper  nulb.  It  is  expected  to  be  useful  in 
many  seismic,  acoustic,  and  electromagnetic  array  processing  applications. 

The  sequence  of  this  section  is  as  follows.  First  we  introduce  the  concept  of  linearly  con¬ 
strained  array  processing.  Then  we  derive  the  optimal  array  weighting  required  to  eliminate  m-1 
interference  sources  incident  on  an  m -element  array  with  a  constrained  look  direction.  The  form 
of  this  optimal  weighting  is  used  as  motivation  for  developing  an  array  processor  based  on  pole- 
zero  filters.  Next,  through  simulations  we  show  how  a  combination  of  both  poles  and  seros  can 
substantially  improve  the  approximation  of  the  optimal  weighting  compared  to  that  obtained 
using  the  conventional  all-zero  method.  Finally,  the  last  part  of  this  section  addresses  the  issues 
of  look-direction  signal  bias  and  inverse  filter  stability  which  arise  because  the  true  output  error  is 
not  being  minimised. 

1.1  Constrained  Array  Processing 

A  constrained  array  processor  such  as  the  one  proposed  by  Frost  [lj,  operates  in  an  environ¬ 
ment  where  a  single  desired  signal  is  incident  upon  an  array  from  a  known  direction  and  several 
interference  signals  are  incident  from  unknown  directions.  The  objective  of  a  Frost  array  proces¬ 
sor  is  to  minimise  the  response  of  the  array  in  the  direction  of  the  interference  signals  while  leav¬ 
ing  the  response  in  the  look  direction  unaltered.  A  block  diagram  of  the  Frost  array  is  shown  in 
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Figure  2.1.  Throughout  this  paper  Z-transform  notation  b  used  to  indicate  a  fixed  linear  digital 
filter.  To  convert  from  the  Z-domain  into  a  function  of  frequency,  the  substitution  z»eJW  b 
made.  The  look-direction  signal  b  assumed  to  be  perpendicular  to  the  array  axis  while  interfer¬ 
ence  signab  are  incident  from  other  unknown  directions.  If  the  look  direction  b  not  perpendicular 
to  the  array  axis,  alignment  filters  usually  in  the  form  of  steering  delays  can  be  added  to  cause 
the  look-direction  signal  to  appear  in  time  coincidence  at  the  output  of  each  delay.  Such  a  sys¬ 
tem  b  referred  to  as  a  signal-aligned  array  [2|.  After  alignment,  the  signab  are  passed  through 
linear  filters  and  then  summed  to  form  the  array  output.  Because  the  look-direction  signal 
appears  identically  at  the  input  to  each  array  filter,  it  will  experience  a  response  through  the 
array  determined  by  the  sum  of  all  the  array  filters.  In  order  to  maintain  an  undbtorted  response 
in  the  look-direction,  thb  sum  b  constrained  to  be  unity.  Subject  to  the  constraint,  the  linear 
filters  are  adjusted  to  minimize  the  output  power  of  the  array.  In  thb  manner,  interference  sig¬ 
nab  are  eliminated  without  dbtorting  the  look-direction  signal.  We  now  determine  the  optimal 
frequency-dependent  weighting  required  to  cancel  m- 1  broadband  directional  interference  sources 
incident  upon  an  m-element  array.  The  analysb  assumes  ideal  propagation  conditions  and  negli¬ 
gible  thermal  noise.  Let  the  Z-transform  of  the  look-direction  signal  be  S(z)  and  the  transform  of 
the  /“  interference  signal  be  A/(z).  The  look-direction  signal  appears  identically  at  all  sensors  of 
the  array  and  b  desired  to  appear  undbtorted  at  the  array  output.  Maintaining  zero  dbtortion  in 
the  look-direction  requires  that 

5(»)  t  "  S{z)  .  (2.1) 

*-i 

Cancelling  S(z)  from  both  sides  of  thb  equation  leaves, 

£  W*(*)  “  I  .  (2-2) 

the  look-direction  constraint.  For  simplicity  we  assume  a  uniformly  spaced  linear  array.  Then 
each  interference  signal  will  experience  a  propagation  delay,  A(l  between  adjacent  sensors  that  b 
determined  by  the  inter-sensor  spacing,  the  speed  of  propagation,  and  the  arrival  angle.  Com- 


plete  cancellation  of  all  interference  signals  requires  that 


*,(*)  E  w*(*)  gAk'1)A‘  "  0  for  l  -  1,2,..., m-1  .  (2.3) 

i 

Equations  (2.1)  and  (2.3)  can  be  combined  and  written  in  matrix  form  as, 


'5(z) 

*i(*)  o 

11...  1 

1  a"A‘  .  .  . 

‘  W,(z)  ' 
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S(z) 

0 

0 

■  i 

1  '  .  .  .  1 

wm(*) 

« 

1 

; 

.  0 

By  induction,  it  can  be  shown  that  the  solution  to  this  set  of  linear  equations  is, 


'wy«) 

>  * 

1 

-£j  all  tingle  delayt 

• 

• 

Use  of  this  frequency-dependent  weighting  in  a  linear  uniformly  spaced  array  wBl  cause  all  of  the 
jammers  to  be  completely  nulled  regardless  of  bandwidth.  However,  since  the  delays  are  usually 
fractional  and  the  frequency  resonances  caused  by  the  poles  are  infinite,  an  exact  rational  realiza¬ 
tion  of  this  weighting  is  not  possible.  Instead,  one  must  be  content  with  a  rational  approximation 
that  is  made  over  a  selected  bandwidth.  Array  performance  using  tapped  delay  lines  operating 
over  a  wide  bandwidth  has  received  some  attention  in  the  literature.  Two  key  papers  are  those 
by  Compton  (3)  and  more  recently  by  Mayhan  and  Simmons  [4|. 

To  gain  insight  into  the  frequency  response  of  the  optimal  weighting,  consider  the  2-element 
single  jammer  case.  When  m— 2,  the  optimal  weight  transfer  functions  are, 
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7igure  2.2  shows  the  frequency  response  of  the  optimal  weights.  The  relevant  frequency  range 
occurs  with  w  between  0  and  it  and  typical  values  of  A  range  between  +  2  and  -2  depending  upon 
the  difference  in  angle  between  the  interference  direction  and  the  look  direction.  The  infinite 
resonance  at  zero  frequency  indicates  that  it  is  physically  impossible  to  maintain  both  unity  d.c. 
gain  in  the  look  direction  and  a  d.c.  null  in  the  interference  direction.  Also  as  the  interference 
direction  approaches  the  look  direction,  A  becomes  small  and  the  frequency  scale  of  the  plot 
expands.  This  causes  the  zero  frequency  resonance  to  dominate  the  0  to  it  frequency  range  and 
makes  approximation  of  the  optimal  weights  more  difficult. 

It  is  informative  to  compare  the  set  of  optimal  weight  equations  given,  in  (2.4)  with  the 
well-known  Wiener  equations.  When  considering  the  effect  of  thermal  noise  on  the  array  output, 
it  turns  out  that  under  certain  circumstances,  the  solution  to  eqs.  (2.4)  will  perfectly  null  out  the 
directional  interference  sources  but  at  the  expense  of  blowing  up  the  thermal  noise.  Thus,  while 
the  signal-to-interference  ratio  is  improved,  there  may  be  a  substantial  degradation  in  the  signal- 
to-thermal  noise  ratio.  For  this  reason,  it  is  important  to  consider  the  effect  of  thermal  noise  in 
order  to  find  a  compromise  between  interference  nulling  and  thermal  noise  suppression.  This  is 
essentially  what  Wiener  theory  does  in  minimizing  the  overall  array  output  power  subject  U,  a 
constraint.  It  can  easily  be  shown  using  Lagrange  multipliers  that  the  weights  which  minimize 
the  mean  square  output  of  the  array  subject  to  a  unity-gain  look-direction  constraint  satisfy  the 
following  equations  [l], 


'l ' 

1 

1 

1 

R„(z)W(z)  * 

• 

/[l  1  .  .  .  lJ-RjJ- 

• 

.1  , 

.1. 

The  mzm  matrix  R„(z)  is  the  spectral  covariance  matrix  of  the  array  signals  after  alignment  and 
W(z)  is  the  vector  of  the  optimal  frequency-dependent  weights. 


look  direction  • 


Figure  2.1  Block  diagram  of  a  Frost  array  with  unity  gain  look-direction  constraint. 


Figure  2.2  Frequency  response  of  optimal  array  weights  for  a  two-element  Frost  array  with  a 
single  interference  source  Incident  on  the  array.  The  Jaimter  arrival  angle  and 
sampling  rate  determine  the  intersensor  delay  i.  (a)  Sensor  1  ( b)  Sensor  2. 


To  connect  this  set  of  equations  to  the  set  given  in  (2.4),  let  (2.4)  be  more  concisely  written  as, 


m 


A(i)W(a) 


Then  multiply  both  sides  of  (8)  by  Ar(z'1)  to  give, 

5(s) 

o’ 

Ar(*-1)A(*)-W(ar)  -  Ar(s‘1)-  '  (2.9) 

I 

o' 

a 

By  interpreting  |S(z)|2  and  |J?|(s)|2  as  the  power  spectra  of  stationary  stochastic  signals,  the 
matrix  Ar(2_1)-A(«)  can  be  interpreted  as  the  spectral  covariance  matrix, 


ArJ,'1)  A(i)  -  R„(s)  . 


(2.10) 


Also,  it  can  easily  be  seen  that, 


AV)  H*(*)l 


(2.11) 


Therefore  eq.  (2.9)  can  be  written  as 


R„(*)-W(z)  -  1 5(«)|! 


Since  eqs.  (2.4)  were  assumed  to  be  satisfied,  it  now  follows  that 


(2.12) 


Z 
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Thos  with  bo  thermal  noise,  eq.  (2.4)  is  equivalent  to  the  Wieaer  sohitioa  given  in  eq.  (2.7). 


2.3  Alternate  Implementation  of  a  Frost  array  processor 

To  simplify  the  adaptivity  issue,  we  introduce  an  alternate  constrained  array  processing 
structure  first  discussed  by  Appiebaum  and  Chapman  (5)  and  more  recently  shown  to  be 
equivalent  to  a  Frost  array  by  Griffiths  and  Jim  [6].  A  block  diagram  of  this  structure  is  shown 
in  Fig.  2.3.  The  bask  idea  is  to  remove  the  look-direction  signal  from  the  adaptive  process  x- 
input.  This  makes  it  impossible  for  the  adaptive  process  to  cancel  signal  components  from  the 
array  output  by  subtraction.  Interference  signals  appear  in  both  the  x  and  d  inputs  and  can 
therefore  be  canceled  from  the  output.  Transformation  between  the  two  array  processors  is 
readily  seen  to  be, 


Vi  (*)' 
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H^*) 
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-1  1  0 

H^z) 

• 

«  -L 
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-1  . 
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• 
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• 
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-1, 

ff-i(*) 

For  the  m«2  case,  this  set  of  equations  becomes. 


Wx{x)  -  V»  -  H(z) 

-  V*  +  U(*)  ■ 


Solving  for  H(z)  gives, 


H(z)  -  l/*  (JW  -  H^z))  , 


and  the  optimal  filter  formulas  of  Eq  (2.6)  reduce  to, 

1  +  r'A 


*(*)  -  V* 


l-r 


(2.14) 


(2.15) 


(2.16) 


(2.17) 


2.4  Simulation  of  a  tw»*«kment  Front  array  processor 

Figure  2.4  compares  the  frequency  response  of  the  converged  Frost  array  filters  with  the  fre¬ 
quency  response  of  the  optimal  weighting.  The  interference  signal  had  a  100%  relative  bandwidth 


L 


J 


Figure  2.3  Alternate  realization  of  e  Frost  errey  using  unconstrained  adaptation. 
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and  an  arrival  angle  of  45  degrees.  The  array  consisted  of  two  elements  and  the  array  filters  each 
had  4  xeros.  The  cross-hatched  dark  lines  indicate  the  optimal  frequency  response  plotted  over 
the  bandwidth  of  the  interference.  The  continuous  curve  represents  the  frequency  response  of  the 
4-xero  Frost  array  filter.  If  the  relative  bandwidth  were  small  enough,  the  optimal  frequency 
response  would  need  to  be  matched  in  magnitude  and  phase  at  only  a  single  frequency.  Such  a 
response  could  easily  be  realized  using  one  complex  weight  as  is  the  case  in  most  narrowband 
adaptive  array  systems.  In  this  paper,  however,  we  are  concerned  with  performance  over  a  wide 
bandwidth.  Figure  2.5  shows  the  broadband  antenna  pattern  formed  by  averaging  narrowband 
patterns  over  the  frequency  range  of  the  interference  signal  The  average  gain  in  the  direction  of 
the  interference  is  20  dB  below  the  signal.  As  more  weights  are  added  to  the  array  filters,  approx* 
imation  of  the  optimal  weight  frequency  response  will  become  more  accurate  and  the  null  depths 
will  become  deeper. 

S.S  A  pole-sero  Front  array  processor 

Aside  from  adding  more  taps  to  the  array  filters,  another  way  to  improve  the  array  perfor¬ 
mance  is  to  use  filters  possessing  both  poles  and  zeros.  In  Fig.  2.2,  we  see  that  the  optimal  weight 
frequency  response  contains  infinite  resonances  at  the  frequencies  n-x/A.  Such  resonances  are 
better  approximated  with  filters  having  both  poles  and  zeros  than  with  zeros  alone.  In  the  next 
section,  we  introduce  a  method  for  adapting  an  array  processor  which  has  array  filters  possessing 
both  poles  and  zeros. 

Development  of  an  adaptive  algorithm  to  adjust  combined  pole-zero  filters  has  been  the  sub¬ 
ject  of  much  research,  but  has  met  with  only  limited  success.  The  reason  for  this  is  that  the 
mean  square  error  function  associated  with  the  pole-zero  filter  is  nonlinear  and  possesses  local 
minima.  This  complicates  the  problem  of  finding  the  global  minimum  of  the  mean  square  error. 
This  problem  does  not  arise  in  an  all-zero  adaptive  filter  since  the  mean  square  error  function  is 
quadratic  in  the  weights  and  the  unique  minimum  can  be  found  by  solving  a  set  of  linear  equa¬ 


tions. 
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Figure  2.4  Frequency  response  of  optimal  array  weights  (cross-hatched  marks) 
compared  to  5  weight  approximation  from  a  Frost  array. 

(a)  Sensor  1  (b)  Sensor  2. 
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Our  approach  bypasses  many  of  the  problems  usually  associated  with  pole-zero  adaptive 
filters.  It  involves  a  reformulation  of  the  error  function  in  a  manner  which  allows  the  polynomial 
associated  with  the  poles  of  the  adaptive  filter  to  be  adapted  in  an  all>zero  form.  This  restores 
the  quadratic  nature  of  the  minimization  problem  to  that  encountered  in  the  all>zero  ease  and 
allows  easy  adaptation  using  Least  Mean  Square  Algorithm  (LMS),  the  Recursive  Least  Squares 
Algorithm  (RLS),  the  sample  matrix  inversion  algorithm  (SMI),  or  any  other  least-squares  algo, 
rithm.  Refer  to  Fig.  2.6  and  assume  for  the  moment  that  the  linear  filters  are  time  invariant. 
Then  the  transfer  function  from  point  d  to  the  output  is  unity  since  the  all-pole  and  all-zero  filters 
cancel.  This  combined  with  the  fact  that  the  look-direction  signal  does  not  appear  in  the  xk 
inputs  allows  the  look-direction  signal  to  pass  through  the  array  undistorted.  The  transfer  func- 

Bklz) 

tion  from  point  xk  to  the  output  is - -  -  and  contains  both  poles  and  zeros.  Instead  of 

1  +  z~lA(z) 

minimizing  the  overall  array  output,  we  minimize  the  signal  appearing  at  the  input  to  the  all- pole 
filter.  In  the  control  literature  this  signal  is  often  referred  to  as  the  equation  error.  By  minimis, 
ing  the  equation  error,  we  are  able  to  adapt  only  all-zero  filters,  yet  effectively  realise  a  pole-zero 
filter  response. 

To  convert  the  system  of  Fig.  2.6  into  a  system  which  uses  constrained  adaptation,  observe 
that  the  transfer  function  from  sensor  k  to  the  output  is  given  by, 


where  £  ni0  —  1  and  the  Nk' (z)  are  unconstrained.  This  leads  to  the  alternate  structural 

kml 

implementation  shown  in  Fig.  2.7.  The  similarity  between  this  system  and  the  Frost  system 
shown  in  Fig.  2.1  is  apparent.  Instead  of  constraining  the  sum  of  all  the  filters  to  be  unity,  only 
the  sum  of  the  first  coefficients  has  been  constrained.  A  look-direction  signal  will  experience  a 


-  19  - 


transfer  function  from  input  to  the  adaptation  error  point  of 

1  +  rl  £  Nk'(z)  (2.20) 

l 

This  response  is  inverted  by  the  all-pole  output  filter  to  generate  the  array  output  and  restore  the 
look-direction  constraint. 

S.fi  Simulation  of  the  pole-iero  Frost  array  processor 

Figure  2.8  compares  the  frequency  response  of  the  optimal  weight  at  each  sensor  with  the 
frequency  response  of  the  converged  2  pole  -  2  zero  filter  obtained  using  the  proposed  pole-zero 
adaptive  array  processor.  As  in  the  previous  simulation,  the  interference  source  had  a  100%  rela¬ 
tive  bandwidth  and  an  arrival  angle  45  degrees  off  the  look  direction.  Notice  the  improvement  of 
this  response  over  the  all-zero  response  plotted  in  Fig.  2.4.  For  further  comparison,  Fig.  2.9 
shows  the  broadband  antenna  pattern  and  the  frequency  response  of  the  array  in  the  signal  and 
interference  directions.  The  interference  null  is  about  40  dB  below  the  look-direction  gain.  This 
represents  a  20  dB  improvement  over  a  conventional  Frost  system  having  the  same  number  of 
degrees  of  freedom.  Fig.  2.10  compares  the  nulling  capabilities  of  the  conventional  Frost  and  the 
pole-zero  Frost  as  a  function  of  jammer  arrival  angle.  Notice  that  as  the  jammer  arrival  angle 
approaches  the  look  direction,  the  null  depth  goes  to  OdB  because  of  the  unity  gain  look-direction 
constraint. 

3.7  Stability  of  the  Inverse  filter 

One  issue  yet  to  be  addressed  is  that  of  stability  of  the  inverse  filter.  In  some  situations,  it 
is  possible  that  the  roots  of  the  polynomial  associated  with  the  inverse  filter  can  move  outside  the 
unit  circle.  When  this  occurs,  the  inverse  filter  will  become  unstable  and  its  output  will  begin  to 
grow  unboundedly.  A  solution  with  poles  outside  the  unit  circle  represents  a  stable  but  noncausal 
filter  and  can  only  be  realized  by  noncausally  processing  the  input  data.  Since  noncausal  process¬ 
ing  is  not  an  available  option  when  operating  in  realtime,  pole  polynomials  with  roots  outside  the 
unit  circle  must  be  avoided. 


The  method  we  propose  for  keeping  all  the  pole  polynomial  roots  inside  the  unit  circle 
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Figure  2.8  Comparison  of  frequency  response  of  optimal  array  Mights  (cross-hatched  marks) 
with  a  2  pole  -  2  zero  approximation,  (a)  Sensor  1  (b)  Sensor  2. 


WIDEBAND 


NORMALIZED  FREQUENCY 


involves  the  addition  of  a  bulk  delay  in  the  desired  response  path  of  Fig.  2.6.  Increasing  the  bulk 
delay  causes  the  roots  of  the  converged  solution  to  lie  closer  to  the  origin.  Choosing  a  delay  too 
large  will  cause  the  resulting  solution  to  lose  quality  and  is  not  desirable.  Choosing  a  delay  too 
small  will  result  in  a  noncausal  solution  and  an  unstable  inverse  filter.  The  idea  is  thus  to  choose 
the  smallest  delay  which  causes  all  the  roots  to  lie  within  the  unit  circle. 

In  the  example  shown  in  Fig.  2.9,  the  converged  solution  was  stable  and  bulk  delay  was  not 
required.  In  fact,  all  of  the  two-element,  single-jammer  experiments  we  simulated  turned  out  not 
to  require  any  bulk  delay.  However,  the  three-element,  two-jammer  cases  tried  required  one  unit 
of  bulk  delay  in  order  to  keep  the  poles  within  the  unit  circle. 

Currently,  we  are  experimenting  with  several  on-line  methods  for  adaptively  choosing  the 
bulk  delay.  A  brute  force  technique  is  to  have  a  series  of  systems  all  operating  on  the  same  input 
data  but  each  using  a  different  delay.  A  simple  criterion  is  used  to  choose  the  stable  system  hav¬ 
ing  the  smallest  delay.  Another  technique  is  to  use  a  single  system  with  a  variable  delay  which  is 
incremented  or  decremented  depending  on  the  stability  of  the  pole  polynomial. 

2.8  Look-direction  signal  bias 

The  Wiener  solution  of  the  conventional  Frost  array  is  not  affected  by  the  presence  of  a 
look-direction  signal.  This  is  readily  apparent  by  referring  to  Fig.  2.3  and  noting  that  the  x-input 
does  not  contain  look-direction  signal.  Therefore  the  look-direction  signal  cannot  contribute  to 
either  the  cross-correlation  between  the  desired-response  input  d  and  the  x-input  or  the  autocorre¬ 
lation  of  the  x-input.  The  Wiener  solution  of  the  pole-zero  Frost  array  is,  however,  affected  by  a 
look-direction  signal.  In  this  section,  we  demonstrate  the  extent  of  this  effect  and  describe  a 
modification  of  the  pole-zero  array  which  eliminates  the  dependence  of  the  Wiener  solution  on  the 
look-direction  signal. 

First,  it  should  be  noted  that  in  situations  where  adaptation  can  be  controlled  so  that  it 
occurs  only  when  the  look-direction  signal  is  not  present,  the  problem  of  the  look-direction  signal 
disturbing  the  Wiener  solution  will  not  exist.  This  is  the  case  in  many  radar  or  spread-spectrum 


applications  where  the  desired  signal  is  pulsed  or  sequenced  in  a  manner  which  is  known  to  the 
receiver.  In  many  continuous  communication  systems  or  seismic  systems,  the  desired  signal  is 
continuous  and  an  alternate  solution  must  be  found. 

To  demonstrate  the  effect  of  the  look-direction  signal  on  the  Wiener  solution  of  the  pole- 
zero  Frost  array,  refer  to  Fig.  2.6  and  assume  that  the  only  source  of  signal  comes  from  the  look 
direction.  In  this  condition,  the  x-input  will  be  zero  and  the  d-input  will  contain  only  signal.  The 
1+  z-1A(z)  filter  will  act  as  a  linear  predictor  and  try  to  reduce  the  look-direction  signal  power. 
If  the  signal  were  narrowband,  the  predictor  would  form  a  notch  filter  at  the  signal  frequency. 
This  transfer  function  would  theoretically  be  inverted  by  the  inverse  output  filter  to  restore  the 
look-direction  signal  at  the  array  output  so  signal  distortion  is  not  the  issue.  The  main  concern  is 
that  degrees  of  freedom  will  be  used  in  attempting  to  internally  null  the  look-direction  signal  and 
thus  cause  a  reduction  of  the  interference  nulling  capabilities  of  the  array.  This  effect  is  referred 
to  as  look-direction  signal  bias. 

One  method  for  alleviating  the  look-direction  signal  bias  problem  is  shown  in  Fig.  2.11.  The 
method  is  based  on  an  idea  proposed  by  Duvall  for  curing  signal  cancellation  in  adaptive  arrays 
[7,8].  For  simplicity,  we  have  shown  only  a  three-element  example.  The  basic  idea  is  to  process 
the  set  of  signals  from  the  sensors  in  a  manner  which  eliminates  the  look-direction  signal  but 
which  keeps  the  phase  relationships  among  the  new  set  of  signals  the  same  as  that  prior  to  pro¬ 
cessing.  The  processed  set  of  signals  is  fed  into  an  adaptive  ’’master”  array  where  the  Wiener 
solution  will  not  be  affected  by  the  look-direction  signal.  The  weights  of  the  master  array  can  be 
copied  into  the  "slave”  array  which  operates  on  the  original  signal  set.  Because  of  the  preserva¬ 
tion  of  phase,  nulls  formed  using  the  slaved  array  will  be  in  the  same  directions  as  those  formed 
in  the  master. 

To  analyze  the  effect  of  preprocessing  on  the  Wiener  solution,  we  define  an  effective  interfer¬ 
ence  environment  that  incorporates  the  preprocessor  into  the  interference  sources.  Each  interfer¬ 
ence  source  will  experience  an  angle-dependent  frequency  transformation  when  passed  through  the 
preprocessor.  The  frequency  response  of  this  transformation  is  determined  by  the  subtraction  of 
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Figure  2.11  Method  For  eliminating  bias  In  Wiener  solution 
due  to  look-direction  signal. 
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signals  appearing  at  adjacent  sensors. 

Again  we  assume  that  the  array  of  sensors  is  linear  and  uniformly  spaced.  A  directional 
interference  source  will  experience  a  delay,  A,  between  two  adjacent  sensors  that  is  determined  by 
the  arrival  angle.  The  frequency  response  of  the  top  subtractor  in  the  preprocessor  is  given  by 

H{u&)  -  1  -  .  (2.21) 

Proceeding  down  the  array,  each  subtractor  will  give  this  same  response  along  with  an  additional 
delay.  The  idea  is  to  use  the  H(u)  transformation  to  modify  the  interference  source  and  concep¬ 
tually  eliminate  the  preprocessor  from  the  array.  As  the  interference  angle  approaches  the  look 
direction,  the  magnitude  of  the  transfer  function  goes  to  zero.  Thus,  the  look-direction  signal  will 
be  removed  and  an  interference  source  appearing  off  the  look  direction  will  be  somewhat 
attenuated  and  high-pass  filtered.  In  this  manner,  the  preprocessor  affects  the  Wiener  solution  by 
modifying  the  weighting  applied  to  the  interference  source. 
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m.  THE  SCALAR  LEAST-SQUARES  LATTICE  FILTER 

3.1  Introduction 

Lattice  filters  were  first  introduced  as  a  means  of  predicting  speech  signals.  Work  on  linear 
prediction  of  speech  using  partial  correlation  techniques  to  whiten  or  decorrelate  the  input  ori¬ 
ginated  in  1969  (9j.  Since  then  the  range  of  applications  for  lattice  filters  has  broadened  to 
include  joint-process  estimation  {10],  signal  enhancement  [llj  and  equalisation  (12).  The  least- 
squares  lattice  filter  that  we  are  now  presenting  was  first  used  by  Satorius  and  Pack  [12]  in  the 
equalization  context.  The  most  interesting  properties  of  the  filter  are  presented  here. 

We  begin  by  presenting  the  orthogonality  principle  of  linear  least-squares  estimation  and  its 
application  to  the  lattice  predictor.  Next,  we  present  the  lattice  predictor  first  used  for  analyzing 
speech  by  decomposing  it  into  a  set  of  orthogonal  components.  We  continue  with  the  structure  of 
the  inverse  of  the  analysis  lattice,  also  known  as  the  synthesis  lattice  in  speech  processing.  We 
also  give  the  necessary  and  sufficient  conditions  for  the  stability  of  that  inverse  filter.  Next,  we 
present  the  joint-process  lattice  for  estimating  a  desired  signal  given  an  auxiliary  one.  Last,  we 
present  simulations  of  the  parallel  modeling  of  an  unknown  all-zero  filter  using  the  joint-process 
lattice. 

3.3  Prediction  and  the  orthogonality  principle 

As  illustrated  in  Fig.  3.1,  the  forward  prediction  problem  is  that  of  estimating  the  next  sig¬ 
nal  value  z(n)  given  M  of  its  past  values,  from  z(n-Af)  to  x(n-l).  Likewise,  backward  predic¬ 
tion  involves  estimating  the  past  value  x(n-M-l)  of  the  signal  given  M  of  its  more  immediate 
past  values  [13]. 

In  the  case  of  least-squares  estimation,  the  orthogonality  principle  [14]  states  that  by  view¬ 
ing  the  signal  samples  as  elements  of  a  vector  space,  the  least-squares  estimate  of  a  sample  is 
found  by  projecting  it  onto  the  subspace  formed  by  the  observed  past  M  samples.  As  a  result, 
not  only  is  the  current  estimation  error  (also  known  as  the  residual  or  innovation)  orthogonal  to 
the  past  M  samples  but  it  is  also  orthogonal  to  all  past  Af  estimation  errors.  Fig.  3.2  illustrates 
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that  orthogonality  principle. 

3.3  The  Lattice  Predictor 

As  shown  in  Fig.  3.3,  the  forward  residuals  Jm* s  are  recursively  formed  by  subtracting  a 
linear  combination  of  the  past  M  samples  of  the  input  from  the  current  sample  z(n).  Likewise 
the  backward  residuals  bm’a  are  formed  by  subtracting  a  linear  combination  of  the  same  M  sam¬ 
ples  from  the  oldest  sample  x(n-M).  The  values  of  the  linear  combination  are  determined  by  the 
reflection  coefficients  hm’a  and  gm’a  to  minimise  the  weighted  sum  of  /£  and  the  weighted  sum  of 

h*. 

Two  properties  are  worth  mentioning  at  this  point.  First,  from  the  orthogonality  principle, 
each  innovation  /  is  orthogonal  to  any  other  innovation  /.  Likewise  the  backward  innovations 
b ’s  are  orthogonal  to  each  other. 

The  second  property  concerns  the  roots  of  the  transfer  functions  from  x  to  Ju  and  from  x 
to  by.  We  can  consider  those  forward  and  backward  prediction  errors  /M(n)  and  bm[n)  to  be  the 
output  of  two  filters  Fm(z)  and  Bm{z)  fed  with  the  same  input  *(n).  Providing  that  km  —  gm, 
(that  occurs  with  stationary  inputs)  the  polynomials  Fm  and  Bm  have  roots  that  are  inverse  of 
each  other  [13]  namely, 

F.M-0  <-> 

Equivalently,  we  can  conclude  that  Fw(?)  has  the  same  coefficients  as  Bm(z)  but  in  reverse 
order. 

3.4  Inverse  of  a  lattice  predictor 

Thus  far  it  seems  that  we  have  not  really  estimated  the  current  sample  i(n)  of  a  signal 
since  we  have  assumed  its  availability  to  construct  the  lattice  predictor.  The  answer  to  that  argu¬ 
ment  has  some  historical  aspects. 

Part  of  the  natural  process  of  speech  generation  is  often  found  to  be  adequately  modeled  by 
an  all-pole  filter  1  /A{z)  fed  with  white  noise  [13].  If  we  use  that  speech  as  the  input  to  a  lattice 


hi  samples  observed  •* 


Fig.  3.1.  Forward  and  backward  predictions  [after  Markel  and 
Gray]. 


Fig.  3.4 

e(n)_L  x(n-l),  x(n-2),  •  •  •  =>  e(n)x  c(n-l).  •  . 

e(n)  is  the  current  estimation  error  or  residual  error.  It  is  also 
known  as  the  innovation  as  it  represents  the  part  of  x(n)  that 
does  not  belong  to  the  subspace  spanned  by  its  past  values,  i.e. 
the  part  that  really  is  new. 
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predictor,  the  forward  innovation  Ju  will  be  white  noise  providing  that  the  order  M  is  sufficient. 
Thus,  the  lattice  transfer  function  from  *  to  fy  is  A(s),  an  estimate  of  the  true  filter  denomina¬ 
tor.  The  process  of  estimating  A(s)  is  known  as  speech  analysis  (13).  Eventually,  we  can  syn¬ 
thesise  part  of  the  speech  by  feeding  white  noise  into  an  inverse  filter  1/A(s)  which  is  an  esti¬ 
mate  of  the  vocal  tract  transfer  function  as  shown  in  Fig.  3.4. 

Thus  the  complete  prediction  process  involves  taking  the  inverse  of  an  all-sero  filter.  We 
now  show  that  the  inverse  can  be  implemented  in  a  structure  that  is  very  similar  to  that  of  a  lat¬ 
tice. 

We  recall  that  the  lattice  recursions  are 

/-(»)  -  /—i(")  +  9m  *— i(»-l)  (3.1) 

*«(»)-*»  /«-i(n)+  *«-i(n-l)  (3.2) 

Since  we  are  interested  only  in  having  /„_ x  in  terms  of  /„,  let  us  rewrite  (3.1)  leaving  (3.2) 
unchanged 

—  / m(n )  -  9m  *«-l(»-l)  (3.3) 

*»(»)  -  *»  fm-1  (»)  +  *»-i(»-l)  (3.2) 

These  recursions  can  be  implemented  as  shown  in  Fig.  3.5  resulting  in  a  structure  having  an 
interesting  similarity  with  the  analysis  filter. 

3.5  Necessary  and  sufficient  stability  conditions  for  the  Inverse 

Although  the  experiment  described  in  Fig.  3.4  clearly  involves  a  stable  l/A(z)  filter,  there 
are  applications  in  which  it  is  important  to  check  the  invertibility  of  the  lattice  predictor.  For 
example  in  pole-zero  filtering  using  the  equation  error,  it  is  necessary  to  invert  the  filter  in  order 
to  restore  the  true  output  error. 

Necettary  conditions 

We  recall  that  the  lattice  recursion  can  be  written  in  terms  of  the  polynomials  F{z )  and 
B{z)  corresponding  to  the  forward  and  backward  time  innovations,  namely  the  forward 


Fig. 3. 3.  Lattice  predictor.  The  /ms  are  the  forward  prediction 
errors  and  the  bm'a  are  the  backward  prediction  errors.  The 
reflection  coefficients  h^'a  and  gm'a  are  updated  to  minimize  both 
the  weighted  sum  of  and  the  weighted  sum  of 
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Fig.  3  4.  Speech  synthesis  process.  The  speech  generation  model  is 
first  identified  by  an  analysis  lattice  predictor.  An  inverse  is  then 
formed  which  yields  the  speech  when  fed  with  white  noise. 


innovation  is  updated  as 


Fm(z)  -  *•„_,(*)  +  gm  i-'B^z)  (3.4) 

Writing  /'m(z)  in  terms  of  its  roots  pm  i  gives 

nd- Pm.i  *'1)  “  +  9m  z'lBm.i(z)  (3.5) 

i  -  0 


The  initial  conditions  for  the  recursion  are  Fo(r)  =»  B^z)  «  1.  Thus  for  m  >  1, 

order  ^Pn.1(s)j  <  m-1 
and 

order  £p„_i(s)j  <  m-1  ■*>  order  J^_15»_x(^)j  <  m 

Thus  the  highest  degree  polynomial  on  the  right-hand  side  of  (3.4)  is  z~l  Bm^(z).  Furthermore, 
the  coefficient  of  the  highest  power  of  *_1  in  s~1Bn.1(s)  is  1.  Equating  the  coefficients  of  the 
highest  degree  terms  of  both  (3.4)  and  (3.5)  gives 

rn-P-.i)-?-  (3-6) 

i  —  l 


Evaluating  the  magnitude  of  both  sides  of  (3.6)  gives 


Thus  if  the  roots  pm  l  have  magnitude  less  than  1,  the  product  of  those  magnitudes  is  less  than  1 
implying  that  the  magnitude  of  the  reflection  coefficient  gm  be  less  than  1  too.  A  necessary  condi¬ 
tion  for  the  roots  of  the  end-polynomial  F^(z)  to  be  within  the  unit  circle  is  therefore 

Jm  !tab"  _>  H<> 


Sufficient  conditions 

Using  a  more  complete  analysis,  Market  and  Gray  [13]  were  able  to  show  that  for  equal 
reflection  coefficients  a  necessary  and  sufficient  condition  for  the  stability  of  1  /Fu(z)  is  that  the 
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reflection  coefficients  for  nil  stages  be  of  magnitude  less  than  1,  namely 

ha  |  <  1,  m  —  1,2,...,  M 

where 

™  ™ 

That  reqnirement  of  equal  reflection  coefficients  is  not  necessarily  stringent  since  in  practice 
an  equivalent  lattice  predictor  meeting  that  condition  is  computed  before  stability  can  be  checked 
(see  section  V). 

3.6  Joint-process  lattice 

The  situation  often  arises  in  signal  processing  when  we  want  to  filter  one  signal  to  identify 
another  signal  or  to  remove  any  correlation  of  it  with  the  former  one.  The  filtering  of  a  to  get  an 
estimate  of  i  can  be  done  by  the  joint-process  lattice  of  Pig.  3.6.  The  backward  innovations  bm 
are  used  as  a  basis  on  which  the  desired  signal  d  is  decomposed.  The  gains  vm's  are  determined 
so  as  to  minimise  the  sum  of  the  t£.  Also  because  of  the  orthogonality,  it  is  possible  to  first  esti¬ 
mate  i  with  v0  bo  yielding  an  error  £q  which  is  then  estimated  by  b ,  yielding  an  error  etc. 

Last,  because  of  the  whitening  that  is  taking  place,  we  expect  the  joint-process  lattice  to 
exhibit  a  convergence  rate  that  is  independent  of  the  z-input  eigenvalues. 

3.7  Simulations 

The  simulation  context  is  the  parallel  modeling  of  an  unknown  all-zero  filter, 
B{z)  mm  l  +  2r_1  +  4*'*  +  2  r-3  +  z~*.  The  joint-process  lattice  is  fed  with  the  same  input  as 
the  unknown  filter  and  the  output  of  the  lattice  is  the  desired  signal  of  the  lattice. 

To  minimise  any  initial  conditions  effect,  we  do  not  look  at  the  convergence  of  the  lattice 
from  its  start-up.  Rather,  we  wait  for  the  weights  to  converge  then  inject  an  impulse  as  a  distur¬ 
bance  at  the  output  of  the  unknown  filter.  Our  criterion  is  therefore  the  time  it  takes  to  reject 
that  disturbance. 


fu(*) 


stable  <  —> 


Furthermore,  to  truly  evaluate  the  convergence  time,  the  simulations  are  noise-free,  yielding 
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a  minimum  mean-square  error  of  the  order  of  -100  dB.  It  is  understood  that  in  actual  applica¬ 
tions  the  noise  level  is  much  higher.  For  evaluation  purposes  however,  only  a  noise-free  environ¬ 
ment  makes  the  algorithm  exhibit  its  full  convergence  potential. 

Fig.  3.7(a)  and  3.7(b)  show  the  mean-square  error  of  the  joint-process  lattice  for  an  input 
eigenvalue  ratio  of  1  (white  noise)  and  10  respectively.  As  expected,  the  convergence  time 
remains  essentially  constant. 


Fig. 3. 7  Disturbance  rejection,  curve  of  scalar  joint-process  lattice 
for  an  input  eigenvalue  ratio  of 
(a)  1 
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IV.  POLE-ZERO  FILTERING  USING  A  TWO-CHANNEL  LATTICE  PREDICTOR 

4.1  Introduction 

In  this  section,  we  present  a  novel  approach  to  adaptive  pole-zero  filtering.  We  show  that 
the  equation-error  method  can  be  implemented  using  a  two-channel  (2-C)  lattice  predictor.  We 
begin  by  presenting  the  2-C  lattice  predictor  as  an  extension  of  the  scalar  predictor  discussed  in 
section  III.  Next,  we  reformulate  the  equation-error  approach  introduced  in  previous  reports  in 
terms  of  linear  prediction  and  then  present  the  equation-error  lattice  predictor  used  for  pole-zero 
filtering.  Last,  a  simulation  involving  the  identification  of  an  unknown  pole-zero  filter  is 
presented  demonstrating  the  relative  insensitivity  of  the  2-C  lattice  predictor  to  an  eigenvalue 
disparity. 

4.2  The  2-C  lattice  predictor 

The  scalar  lattice  predictor  presented  in  section  m  can  be  generalized  to  a  multi¬ 
dimensional  lattice  predictor  when  the  process  to  be  whitened  is  itself  multi-dimensional.  The 
algorithm  is  detailed  in  [15|.  For  the  moment,  we  restrict  ourselves  to  considering  only  the  2-C 
case.  The  multi-dimensional  case  is  Considered  in  section  V  for  multi-element  antenna  array  pro¬ 
cessing. 

Given  a  vector 


the  lattice  predictor  constructs  a  sequence  of  forward  and  backward  innovation  vectors 


In  particular,  fu  (see  Fig.  4.1)  is  the  forward  prediction  error  of  Xj  gt  the  past  M  values  of  xx 
and  x2,  namely, 

fit  “  *i(")  ~  *i(")  |  . *i(n-M),  x2(n-l),...,  x2(n-Af) 

Recognizing  the  forward  prediction  error  as  such  will  allow  us  to  use  the  2-C  lattice  predictor  in 


pole-zero  array  processing. 


4.3  Reformulation  of  the  equation  error  approach 

The  equation  error  method  has  been  discussed  in  detail  in  reference  [16],  The  objective  is  to 
filter  the  auxiliary  input  z  with  an  autoregressive  and  moving  average  (ARMA)  filter, 
B(z)/(  1  +  z~lA(z)),  i.e.  using  poles  and  zeros  so  that  the  output  error  is  minimized  (see  Fig. 
4.2(a)).  As  shown  in  the  previous  reports,  the  A(z)  and  B(z)  filters  can  equivalently  estimated  by 
minimizing  the  equation  error  between  d  filtered  through  (1  +  z~lA(z))  and  z  filtered  throuf*' 
B(z)  (see  Fig.  4.2(b)). 

As  shown  in  Fig.  4.2(c),  the  equation  error  method  can  be  thought  of  as  the  minimization  of 
the  forward  prediction  error  of  d  given  the  past  values  of  d,  the  current  value  of  z  and  the  past 
values  of  z.  This  is  essentially  what  a  2-C  lattice  predictor  does. 

4.4  Equation-error  lattice  predictor 

We  are  changing  the  seemingly  scalar  problem  of  pole-zero  filtering  into  a  multivariable 
problem  by  embedding  the  primary  and  auxiliary  signals  into  a  vector. 

We  recall  that  in  the  2-C  lattice  predictor,  the  prediction  of  any  element  of  the  input  vector 
[z,  x2j  is  based  upon  past  values  of  both  elements  zl  or  z2  of  that  vector.  The  equation  error 

approach  requires  that  we  also  base  the  prediction  of  d  on  the  current  value  of  z.  Therefore  the 
d-input  needs  to  be  delayed  by  one  in  order  to  obtain  a  correct  prediction  as  shown  in  Fig.  4.3. 

As  in  the  scalar  case,  we  expect  this  novel  adaptive  scheme  to  be  relatively  insensitive  to 
any  disparity  in  the  eigenvalues  of  the  input  vector.  Such  a  property  might  be  very  useful  in 
some  instances.  In  actual  system  identification  for  example,  the  auxiliary  signal  z  (the  input  to 
the  unknown  plant)  is  seldom  white  thus  has  an  eigenvalue  spread.  But  even  if  it  is,  the  primary 
signal  d  is  nothing  but  colored  since  it  is  the  plant  output. 


1% 


VS2l 


Fig.  4  5  Equation  error  approach  using  a  2-C  lattice  predictor  of 
order  M. 

The  output  of  the  predictor  is 

ft(T i)  =  d(n)-d(n)  d(n-l),:..,  d(n-Jlf),  z(n ) . x(n-M  + 1) 

The  transfer  function  from  1  to  4  and  from  3  to  4  identify  the  poly¬ 
nomials  A(z)  and  B(z)  namely, 

#n(z)  =  z_I(l  +  z_li4(z))  and  //^(z)  =  z~lB{z). 


4.5  Simulations 


The  simulation  context  is  the  modeling  of  an  unknown  4th  order  Butterworth  filter,  a  filter 
with  both  poles  and  zeros.  The  predictor  z-input  is  the  input  to  the  Butterworth  filter.  The 
predictor  i-input  is  the  output  of  the  Butterworth  filter. 

As  in  section  m,  we  only  look  at  the  time  to  reject  an  impulse  disturbance  at  the  output  of 
the  unknown  filter.  By  doing  this,  we  minimize  the  effect  of  initial  conditions.  Also,  the  simula¬ 
tion  is  noise-free  in  order  to  evaluate  the  full  convergence  potential  of  the  algorithm. 

Fig.  4.4(a)  and  4.4(b)  show  the  mean  square  prediction  error  of  the  2-C  lattice  for  an  z- 
input  eigenvalue  ratio  of  1  (white  noise)  and  10  respectively. 

The  convergence  time  is  slightly  longer  for  a  ratio  of  10  than  for  pure  white  noise  (650  itera¬ 
tions  vs.  550)  but  it  remains  very  rapid. 


V.  TWO-ELEMENT  ANTENNA  PROCESSING  USING  A  POLE-ZERO  LATTICE 


5.1  Introduction 

The  equation  error  approach  as  applied  to  antenna  array  processing,  requires  the  formation 
of  the  inverse  filter  1  /[I  +  z’Mfz))  in  addition  to  the  estimation  of  A(z)  and  B(z).  In  this  sec¬ 
tion,  we  first  present  a  way  of  forming  that  inverse  and  checking  its  stability.  Next,  we  show  an 
implementation  of  the  pole-zero  lattice  predictor  in  the  adaptive  array  configuration  presented  in 
section  II  of  this  annual  report.  By  removing  the  constraints  in  the  antenna  adaptive  filters,  the 
equation  error  array  processing  structure  is  suitable  for  implementation  using  the  rapidly  converg¬ 
ing  lattice  predictor  we  are  considering.  Last,  we  present  the  simulation  results  for  a  wideband 
jammer  cancellation  experiment. 

5.2  Inverse  of  1  -I-  r1jt(i) 

To  restore  the  true  output  error,  we  need  to  be  able  to  invert  the  1  +  z~l  A(z)  polynomial, 
i.e.  the  transfer  function  from  2  to  4  in  Fig.  4.1.  This  requires  that  H14[z)  mm  1  +  z~l  A(z)  has  a 
stable  inverse.  But  we  need  to  invert  only  that  transfer  function  unlike  in  section  II  where  we 
inverted  the  complete  upper  path  of  the  forward  innovations.  Doing  this  also  inverts  H34(z)  and 
thus  requires  the  latter  to  be  minimum-phase.  Such  a  double  constraint  is  certainly  sufficient  to 
ensure  the  stability  of  1  /(I  +  z"lA(z))  but  is  not  necessary. 

A  computationally  more  complex  method  for  checking  the  stability  of  the  polynomial  Hu(z) 
is  to  convert  it  to  a  scalar  lattice  form  using  what  is  known  as  the  ’’step-down"  procedure  [13|. 
Then  we  need  only  recall  that  the  necessary  and  sufficient  conditions  for  Hu(z)  to  have  roots 
within  the  unit  circle  is  that  all  reflection  coefficients  of  its  lattice  implementation  be  of  magni¬ 
tude  less  than  one. 

6.3  Pole-sero  adaptive  lattice  for  array  processing 

The  array  processor  shown  in  Fig.  5.1  is  a  pole-zero  extension  of  the  Griffiths  and  Jim  gen¬ 
eralized  sidelobe  canceller  [<J|.  It  is  presented  in  detail  in  section  2.3  of  this  report.  Any  signal 


coming  from  the  look  direction  ia  constrained  to  pass  through  the  array  undistorted.  The  prepro¬ 
cessing  of  the  sensor  outputs  allows  the  constraints  to  be  removed  from  the  adaptive  algorithm 
and  allows  us  to  use  the  same  lattice  predictor  as  we  used  in  the  system  identification  experiment. 

As  suggested  previously,  a  scalar  lattice  predictor  equivalent  to  1  +  s~1A(z)  needs  to  be 
computed  and  the  magnitude  of  the  reflection  coefficients  checked  before  an  inverse  is  actually 
implemented. 

5.4  Simulations 

Fig.  5.2  shows  the  simulation  results  of  the  pole-zero  adaptive  lattice  used  for  a  2-element 
antenna.  Those  results  are  to  be  compared  with  the  simulations  presented  in  section  2.6  of  this 


Fig  5.1  Array  processing  using  a  pole-zero  adaptive  lattice. 

Note  that  the  output  is  delayed  by  one  time  unit  due  to  the  delay 
in  the  primary  input  path. 
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Fig.S.z  2-C  lattice  antenna  array. 

(a)  Broadband  antenna  sensitivity  pattern. 

(b)  Frequency  response  in  the  jammer  direction. 

(c)  Frequency  response  in  the  signal  direction. 
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VI.  MULTI-ELEMENT  ARRAY  PROCESSING  USING  A  POLE-ZERO  LATTICE 

6.1  Introduction 

The  use  of  the  equation  error  approach  in  multi-element  array  processing  was  addressed  in 
section  II  of  this  report.  It  is  summarized  here  for  convenience  and  in  order  to  demonstrate  the 
connection  to  dimensional  lattice  prediction. 

As  shown  in  Fig.  6.1(a),  the  pole-zero  filters  Br(z)l(  1  +  z'1  A,  (z))  need  to  be  adjusted  to 
minimize  the  output  error  which  is  similar  to  minimizing  the  equation  error  shown  in  Fig.6.1(b). 
Equation  error  minimization  is  equivalent  to  minimizing  the  prediction  error  of  d  given  past 
values  of  i,  the  current  values  of  z/s  and  the  past  values  of  x/s  (see  Fig.  6.1(c)). 

6.2  Implementation  of  the  multi-channel  pole- two  lattice 

Similarly  to  what  was  done  in  the  2-C  ease,  multi-channel  pcle-zero  filtering  can  be  embed¬ 
ded  into  multi-channel  prediction.  The  implementation  of  such  a  predictor  in  lattice  form  is 
shown  in  Fig.  6.2.  Again,  we  require  that  the  primary  input  path  have  a  unit-delay  in  order  to 
base  the  estimation  of  the  primary  signal  in  the  current  values  of  the  auxiliary  signals  as  well. 

6.3  Limitations  of  the  multi-channel  pole-aero  lattice 

Although  theoretically  appealing,  the  multi-dimensional  pole-zero  lattice  requires  a  large 
number  of  computations.  In  particular,  there  are  2  [(P  +  1)  X  (P  +  1)]  matrix  inversions  each 
time  a  stage  is  updated  where  P  +  1  is  the  total  number  of  lattice  inputs. 

The  forward  and  backward  innovations  are  updated  from  order  m-1  to  order  m  as  follows 
[15| 

/«(«)- /«.-i(«)+  <7«.(n-l)*«,-i(»-l) 

Mn)  —  f/«,(n-l)/*_1(n)  +  bm.i(n- 1) 
where  the  crossover  gain  matrices  Gm  and  Hm  are  computed  by 

<7*(n-l)  -  /tf(»-l)££M«-2) 


(1  +  *-%(«)) 


1  + 


Bp{t) 


1  ♦  t"Mp(*) 


Fig  4.1  Equation  error  in  terms  of  prediction  error 

(a)  Multi-element  array  Altering  with  poles  and  zeros  to  minimize 
the  output  error  e. 

(b)  Processing  using  only  2eros  to  equivalently  minimize  the  equa¬ 
tion  error  s'. 

(c)  e'  is  in  fact  the  prediction  error  o£  d  given  past  d's,  the 
current  xp's  and  the  past  xp's.  The  filter  A  (z)  is  defined  as 

A[e)  =  -1  +  f[  (1  +  z"%(z)) 
p  *  i 

The  filters  Bp{z)  are  defined  as 


=  ^{z)  ft  (1  +  z'%(*)) 


Fig.  6.2. .  Equation  error  approach  using  a  multi-C  lattice  predictor. 
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where  the  ((P  +  1)  X  (P  +  l)j  matrices  E^!x  and  E^lx  are  the  inverses  of  the  minimized  sums 
of  the  backward  and  forward  prediction  errors  respectively.  Those  inversions  take  on  the  order  of 
P 2  operations  which  make  the  crossover  gain  computations  on  the  order  of  P3  operations.  The 
overall  lattice  recursion  is  then  dominated  by  a  count  on  the  order  of  M  X  P3  operations  where 
M  is  the  number  of  stages.  If  P  is  small  with  respect  to  M  that  computational  disadvantage  may 
not  be  too  crucial  but  if  P  is  large  or  simply  close  to  M,  any  improvement  in  convergence  speed 
brought  by  the  lattice  would  be  heavily  paid  in  algorithmic  complexity. 
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vn.  CONCLUSIONS 

In  this  report,  we  have  presented  a  new  wideband  adaptive  array  processing  structure 
Using  a  derivation  of  the  optimal  wideband  array  weighting  function,  we  conjectured  that  filters 
possessing  both  poles  and  zeros  should  be  capable  of  improving  the  performance  of  the  adaptive 
array.  Using  the  "equation-error”  approach,  we  developed  an  adaptive  method  for  adjusting  the 
poles  and  zeros  of  the  array  filters.  Through  simulations,  we  showed  that  the  equation  error 
approach  indeed  does  improve  the  nulling  capabilities  of  the  array. 

Also  in  this  report,  we  discussed  scalar  and  multi-channel  lattice  filters.  Applications  of  the 
lattice  to  linear  prediction,  joint-process  estimation,  general  equation-error  pole-zero  filtering,  and 
pole-zero  array  processing  were  discussed.  The  major  motivation  for  the  use  of  lattice  filters  is 
their  fast  convergence  properties.  Simulations  were  presented  which  confirmed  this  fact. 

The  information  presented  in  this  report  is  by  no  means  the  final  word  on  equation-error 
array  processing.  Several  issues  need  further  consideration  a  few  of  which  w»  are  listed  below. 
Maintaining  stability  of  the  inverse  output  filter  is  a  major  issue.  Techniques  to  easily  check  this 
filter  possibly  using  the  lattice  approach  need  further  investigation.  A  technique  for  adaptively 
choosing  the  minimum  bulk  delay  which  stabilizes  the  pole  polynomial  needs  to  be  further 
developed  and  analyzed.  Additional  experiments  should  be  conducted  with  the  method  presented 
for  elimination  of  look-direction  signal  bias.  The  algorithm  for  implementing  a  multi-channel  lat¬ 
tice  filter  should  be  modified  in  order  to  reduce  the  computational  complexity  from  being  propor¬ 
tional  to  the  cube  of  the  number  of  channels  down  to  the  square  of  the  number  of  channels. 
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